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COFDM DEMODULATOR WITH 
AN OPTIMAL FFT ANALYSIS WINDOW POSITIONING 

The present invention relates to a Coded Orthogonal Frequency Division 
5 Multiplex (COFDM) demodulator. 

Fig. 1 is intended to illustrate the principle of a COFDM modulation. Data 
packets to be transmitted are put in the form of N complex coefficients associated with N 
respective frequencies (or carriers). Number N of the frequencies is equal, for example, 
to 1705 for the so-called "2K" mode and to 6817 for the so-called "8K" mode, in digital 
10 television radio transmission. Each complex coefficient corresponds to a vector which is 
illustrated in Fig. 1 as starting from a frequency axis at a point indicating the frequency 
associated with the coefficient. 

The set of these N coefficients is processed by inverse fast Fourier transform 
(IFFT), which generates a "symbol" formed of a sum of modulated carriers, each carrier 
15 having an amplitude and a phase determined by the associated complex coefficient. The 
symbol thus generated is transmitted. 

Conventionally, in radio transmission, the width of the transmission channel is 6, 
7, or 8 MHz and each carrier is separated from the next one by a frequency difference Af 
= 1/Tu. Tu is the transmit time of a symbol and is called the useful duration. The useful 
20 duration, is on the order of 224 [is in 2K mode and 896 |is in 8K mode, for a 8-MHz 
passband. 

Upon reception, a receiver submits the symbol to the inverse processing, that is, 
mainly, a fast Fourier transform (FFT) to restore the initial complex coefficients. 

As shown in Fig. 1, certain regularly distributed vectors PI, P2, P3... have a 
25 known constant value. These vectors, or the corresponding carriers, are said to be pilot. 
They are used to reflect the distortions undergone by the transmitted channel and by the 
information that they provide on the channel response, they enable correcting the 
unknown vectors located between the pilots. 

Fig. 2 illustrates a transmission of several successive symbols Sn-1, Sn... As 
30 shown, each of these symbols is preceded by a guard interval Tgn-1, Tgn, which is a 
copy of a portion of the end of the corresponding symbol. The guard intervals are often 
defined by a fraction of useful time period Tu. Conventional values of the guard interval 
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are Tu/32, Tu/16, Tu/8, or Tu/4. 

The guard intervals are used to avoid inter-symbol modulation distortions caused 
by an echo of the transmission. Fig. 2 also shows an echo SEn-1, TgEn-1... of the 
transmitted signal. This echo is delayed with respect to the main signal by a time period 
5 shorter than that of a guard interval Tg. 

Each symbol S is normally analyzed by the FFT circuit of the receiver in a 
window W, or FFT analysis window, of same length as the symbol. If there was no 
guard interval, an analysis window W coinciding with a symbol of the main signal would 
be astride two symbols of the echo signal. This would cause an error which would be 
10 difficult to correct in the calculation of the FFT Fourier transform. 

Guard interval Tg, provided that it is greater than the delay or the advance of the 
echo, provides an adjustment margin for analysis window W so that it only coincides 
with portions belonging to the same symbol, in the main signal as well as in the echo. 
The fact for an analysis window to be astride a symbol and its guard interval introduces a 
15 phase shift which is corrected by means of the above-mentioned pilots. 

In Fig. 2, symbol Sn-1 must be analyzed in a window Wn-1 of duration Tu that 
can be placed indifferently in a window delimited by times ta and tb, time ta 
corresponding to the beginning of the guard interval of echo TgEn-1, and time tb 
corresponding to the end of symbol Sn-1. Similarly, symbol Sn must be analyzed in a 
20 window Wn of duration Tu that can be positioned indifferently in the window delimited 
by times tc and td, time tc corresponding to the beginning of the guard interval of echo 
TgEn, and time td corresponding to the end of symbol Sn. 

Fig. 3 schematically shows the place of pilots in the symbols. The symbols are 
gathered in frames of 68 symbols, conventionally in digital television radio transmission 
25 (standard ETSI EN 300 744, VI. 4.1). 

In Fig. 3, each line represents a symbol and each box represents the position of a 
carrier. The carriers are defined as going from a position 0 to a position Kmax, Kmax 
being equal to 1704 in 2K mode and 6816 in 8K mode. Indeed, only a portion of the 
possible frequencies is used, especially due to risks of losses at the channel border. The 
30 pilots are of two types. 

On the one hand, there are, in each symbol, continuous pilots Pc. The continuous 
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pilots correspond to specific frequencies distributed in the channel. In the above- 
mentioned ETSI standard, there are 45 in 2K mode and 177 in 8K mode. The continuous 
pilots are present in all the symbols and always take up the same frequency position. In 
Fig. 3, only the continuous pilots corresponding to positions 0 and Kmax have been 
5 shown. 

On the other hand, there are, in each symbol, so-called "scattered pilots" Pr, 
which are arranged every 12 carriers, and shifted by three positions between two 
successive symbols. Thus, every four symbols, the same arrangement of scattered pilots 
Pr can be found. 

10 Initially, the FFT analysis windows are roughly positioned, for example, by a 

method of intercorrelation of the received signal. The continuous and scattered pilots, of 
constant amplitude on transmission, are then used to finely position analysis windows W. 

For this purpose, at the receiver, the complex time received signal, after having 
been put in baseband, is provided to a fast Fourier transform unit providing the symbol in 

15 the frequency field. The pilots are sampled from this symbol. The sampled pilots enable 
estimation of the frequency response of the channel which, after having undergone an 
inverse Fourier transform, provides the estimation of the channel pulse response. The 
estimation of the channel pulse response is used to finely position analysis window W. 

However, a problem exists due to the fact that the pilots provide but an 

20 incomplete description of the channel. Indeed, a taking into account of several 
successive symbols (at least 4) only enables having an image of the channel every three 
points. As a result, the estimation of the pulse response of the channel obtained from the 
pilots exhibits a periodization of period equal to useful duration Tu divided by three. 

Fig. 4 A shows an example of a curve 1 representative of the module estimation of 

25 the frequency response of a transmit channel. Black disks 2 correspond to the points of 
curve 1 available from the continuous and scattered pilots of four successive symbols. 
The estimation of the pulse response of the transmit channel is obtained by the inverse 
Fourier transform (IFFT) of the points of the estimation of the frequency response 
obtained for the continuous and scattered pilots. 

30 Fig. 4B shows the module estimation of the pulse response of the channel 

obtained from the points of the module estimation of the frequency response shown by 
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the black disks in Fig. 4A. The module estimation of the pulse response is schematically 
shown by a periodic series of pulses 4 having a period equal to useful duration Tu 
divided by 3. For each period, a pulse corresponds to the main path taken by the signal 
and the other pulses correspond to echoes. In practice, the estimation of he pulse 
5 response corresponds more to a periodic series of "blunted" peaks. 

Fig. 5 schematically illustrated the steps of a conventional method for positioning 
an analysis window of the FFT circuit based on the module estimation of the pulse 
response. The method comprises the steps of: 

searching, over a period of the module estimation of the pulse response, pulse 7 
10 of maximum amplitude which corresponds to the main path; 

displacing a window FE, having its width corresponding to guard interval Tg, 
with respect to main path 7 from an initial position 8 (shown in dotted lines) to a final 
position 9 (shown in full line) and determining, for each position of window FE, the 
"energy" of the module estimation of the pulse response in window FE, the "energy" for 
15 example corresponding to the sum of the amplitudes present in window FE; and 

refining the positioning of analysis window W from the position of window FE 
corresponding to the maximum "energy". 

The position of window FE corresponding to a maximum energy is generally 
further used to determine the estimation of the frequency response for the carriers other 
20 than continuous or distributed carriers, represented by white disks 6 in Fig. 4A. For this 
purpose, an interpolation filter is generally applied to the estimation of the pulse 
response, the positioning of which is refined according to the maximum energy position 
of window FE. 

Such methods to refine the positioning of the FFT analysis window and the 
25 positioning of the interpolation filter are implemented on the demodulator referred to as 
STV0360 sold by the applicant. 

In the previously-described method for positioning the FFT analysis window, the 
width of window FE is equal to guard interval Tg. This means that all the echoes outside 
the guard interval are not taken into account to refine the positioning of the FFT analysis 
30 window. However, echoes are generally present outside of the guard interval. The fact 
for them not to be taken into account to refine the positioning of the FFT analysis 
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window may alter the quality of the demodulation. However, the width of window FE 
cannot be increased beyond the guard interval since this risks making the FFT analysis 
window positioning method unstable. 

An object of the present invention is to provide an optimal positioning method of 
the FFT analysis window for a COFDM demodulator enabling taking into account 
echoes outside the guard interval. 

To achieve this object, the present invention provides a method of COFDM 
demodulation of a signal received from a transmit channel, comprising the steps of 
performing the fast Fourier transform of the received signal in a window corresponding 
to a symbol, each symbol comprising several carriers modulated in phase and/or in 
amplitude, some of which are pilots, and being adjacent to a guard interval reproducing 
part of the symbol; providing a set of estimated values of the pulse response in module 
from the pilots; determining coefficients, each coefficient being obtained from the 
product of said set and of a filtering function for a determined relative position of the 
filtering function with respect to said set; determining the maximum coefficient and the 
corresponding relative position; and positioning said window according to said relative 
position corresponding to the maximum coefficient, the filtering function comprising a 
central portion of constant amplitude and of duration equal to the duration of the guard 
interval, surrounded with non-zero decreasing sides. 

According to an embodiment of the present invention, the step of providing the 
set of estimated values of the pulse response in module comprises a step of provision of a 
set of estimated values of the frequency response of the transmit channel based on the 
pilots, and a step of transformation of said set of estimated values of the frequency 
response by inverse fast Fourier transform. 

According to an embodiment of the present invention, for each relative position 
of the filtering function with respect to the set of estimated values of the pulse response 
in module, the coefficient is determined based on the sum of the products of the 
estimated values of the pulse response in module and of the filtering function. 

According to an embodiment of the present invention, the set of estimated values 
of the pulse response in module is periodic, the filtering function having a total width 
smaller than the period of the set of estimated values of the pulse response in module. 
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According to an embodiment of the present invention, the filtering function is a 
stepped function. 

According to an embodiment of the present invention, each side comprises at 
least a first and a second adjacent steps each having a non-zero amplitude strictly smaller 
5 than the amplitude of the central portion, the double of the amplitude of the second step 
being greater than the sum of the amplitude of the central portion and of the amplitude of 
the first step. 

According to an embodiment of the present invention, the duration of the second 
step is smaller than the duration of the first step. 
10 According to an embodiment of the present invention, the filtering function is 

symmetrical. 

The present invention also provides a COFDM demodulator intended to receive a 
signal received from a transmit channel, comprising a circuit of fast Fourier transform of 
the received signal in a window corresponding to a symbol, each symbol comprising 

15 several carriers modulated in phase and/or in amplitude, some of which are pilots, and 
being adjacent to a guard interval reproducing part of the symbol; a circuit for providing 
a set of estimated values of the pulse response in module from the pilots; a circuit for 
determining coefficients, each coefficient being obtained from the product of said set and 
of a filtering function for a determined relative position of the filtering function with 

20 respect to said set; a circuit for determining the maximum coefficient and the 
corresponding relative position; and a circuit for positioning said window according to 
said relative position corresponding to the maximum coefficient, the filtering function 
comprising a central portion of constant amplitude and of duration equal to the duration 
of the guard interval, surrounded with non-zero decreasing sides. 

25 The foregoing object, features, and advantages, as well as others, of the present 

invention will be discussed in detail in the following non-limiting description of specific 
embodiments in connection with the accompanying drawings. 

Fig. 1 , previously described, shows phase- and amplitude-modulated carriers in a 
COFDM transmission system; 

30 Fig. 2, previously described, shows signals received by a COFDM demodulator 

and windows of FFT analysis of the signals; 
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Fig. 3, previously described, schematically shows the position of pilots in 
symbols; 

Figs. 4A and 4B, previously described, schematically show examples of the 
module estimations of the frequency response and of the pulse response of the transmit 
5 channel; 

Fig. 5, previously described, illustrates an example of a conventional method for 
positioning the FFT analysis window; 

Fig. 6 shows an example of the forming of a demodulator according to the present 
invention; 

10 Fig. 7 shows an example of the forming of a filtering function implemented by 

the method according to the present invention of FFT analysis window positioning; 

Figs. 8 and 9 illustrate steps of the method according to the present invention for 
positioning the FFT analysis window; and 

Figs. 10A, 10B to 13 A, 13B illustrate certain advantages of the method according 
15 to the present invention of FFT analysis window positioning. 

Fig. 6 shows an example of a COFDM demodulator according to the present 
invention. The received signal comprises continuous pilots, scattered pilots, and data 
carriers. 

In Fig. 6, an input E of the demodulator receives a signal IF of intermediary 
20 frequency enabling sampling, for example, 36 MHz. Signal IF corresponds to the signal 
received after various frequency changes or transpositions. 

Input E is coupled to an analog-to-digital converter 10 (ADC) which digitizes 
input signal IF. Analog-to-digital converter 10 drives a frequency change unit 12. Unit 
12 provides a signal substantially in baseband, the signal spectrum at the output of unit 
25 12 being centered on a frequency substantially equal to zero. Unit 12 is coupled to a unit 
14, enabling on the one hand fine setting of the central frequency of the signal spectrum 
and, on the other hand, providing time samples at times appropriate to the subsequent 
processing. At the output of unit 14, the signal spectrum is centered on a frequency equal 
to 0 and the number and the time position of the samples are adapted to the 
30 transformation by Fourier transform which is performed in the next unit. Unit 14 is 
controlled by connections 15 and 15 ! connecting unit 14 to a unit 16 for processing the 
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continuous and scattered pilots. 

The output of unit 14 drives a fast Fourier transform unit 20 (FFT) which 
provides the frequencies corresponding to a symbol. Unit 20 is driven by a unit 22 which 
provides, via a connection 24, a signal for setting the analysis window of the Fourier 
5 transform. 

The output of unit 20 is coupled to unit 16 which performs the extraction and the 
processing of the continuous and scattered pilots. Unit 16 provides on connections 15 
and 15' the signals intended to correct the central frequency of the spectrum and of the 
sampling frequency of the signal. 

10 The output of unit 20 drives a unit 30 in which the signal is connected by means 

of an estimation of the frequency response of the signal. The estimation of the frequency 
response of the channel is performed in unit 16 by means of the pilots. This estimation is 
provided by unit 16 on a connection 35, having a branch 35a coupled to unit 30. At the 
output of unit 30, the signal especially comprises the carriers conveying the data. 

15 The estimation of the frequency response of the channel, provided by unit 16, 

supplies, via connection 35 and a branch 35b of connection 35, an inverse fast Fourier 
transform unit 36 (EFFT), to determine the pulse response of the channel. Unit 36 
provides the pulse response of the channel to unit 22, to dynamically adjust the 
positioning of the FFT analysis window. Unit 22 is connected to unit 30 to dynamically 

20 adjust the position of an interpolation filter used to determine the estimation of the 
channel frequency response for carriers other than the continuous and scattered pilots. 

The processing of the carriers conveying the data is ensured in a data processing 
and provision circuit 40. Circuit 40 has a conventional structure and may comprise, as 
shown in Fig. 4, a symbol disinterlacing unit 42, a so-called "demapping" unit 44, a bit 

25 disinterlacing unit 46, and an error correction unit 48 (FEC). The output of unit 48 forms 
output S of circuit 40 and of the demodulator and provides data corresponding to the 
transmitted data. 

The method to refine the positioning of the FFT analysis window comprises the 
previously-described steps. The present invention consists of using a specific filtering 
30 function instead of previously-described window FE enabling taking into account certain 
echoes outside of the guard interval. 
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Fig. 7 shows an example of a filtering function FE according to the present 
invention. The scales, especially according to the time axis, are not respected. 

In the present example of implementation, filtering function FE comprises N 
successive steps. Each step of index j, j=l to N, is characterized by a duration AT; and 



pulse response with respect to filtering function FE may be performed from a portion of 
the module estimation of the pulse response of the channel of duration equal to the period 
of the module estimation of the pulse response, that is, equal to useful duration Tu 
divided by 3 in the present example. The total sum of durations ATj, j=l to N, can then 

10 be greater than the period of the module estimation of the pulse response. However, the 
determination of the energy of the module estimation of the pulse response with respect 
to filtering function FE can be performed based on the entire module estimation of the 
pulse response. The total sum of durations ATj, j=l to N, is then smaller than the period 

of the module estimation of the pulse response. 
15 Amplitudes Lj increase from amplitude Lj of the first step to a maximum 

amplitude LMAX corresponding to a step of index mi, then decreases from maximum 
amplitude LMAX to amplitude L^ of the last step of index N. Duration AT m [ of the 

step of index mi is equal to the duration of guard interval Tg. The filtering function can 
be considered as having a constant central portion, the step of index mi, surrounded with 
20 decreasing sides, the steps of indexes j=l to mi-1 and j=mi+l to N. The decreasing sides 
are non-zero since at least one of amplitudes Lj for j=l to mi-1 and for j=mi+l to N is 

different from zero. 



5 by an amplitude Lj. The determination of the energy of the module estimation of the 



Durations ATj, j=l to N, are such that: 




25 



This means that the step of index mi is substantially located at the center of 



filtering function FE. 

Amplitude Lj of the step of index j follows the following relation, for j=l to N: 

2xL j >Lj. 1 +Lj +1 
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As an example, the following relations apply: 
ATj+i > ATj at least for mi < j < mi+4 

ATj+i < ATj at least for mi-4 < j < mi 

More specifically, advantageously, the following relations apply: 
5 ATj+i > 1 .5 x ATj at least for mi < j < mi+4 

ATj+i < 1.5 x ATj at least for mi-4 < j < mi 

Advantageously, number N of steps is odd. As an example, 19 steps are shown in 
Fig. 7. Further, filtering function FE is advantageously symmetrical, which translates, 
for j=l to N, as the following relations: 
10 ATj=AT N +!.j 

Lj=L N +i_j 

Generally, for each index, amplitude Lj and duration ATj are set according to the 

used modulation, especially according to the duration of guard interval Tg. Further, 
whatever the used modulation, amplitudes Li and L>j respectively associated with the 

15 steps of indexes 1 and N are equal. 

Figs. 8 and 9 schematically illustrate two steps of the method for positioning the 
FFT analysis window according to the present invention implementing filtering function 
FE shown in Fig. 7. The scales, according to the time axis, are not respected. As 
described previously, to determine the energy of the estimation of the pulse response 

20 with respect to filtering function FE, a portion only of the module estimation of the pulse 
response having a duration equal to a period of the module estimation of the pulse 
response. Said portion of the estimation of the pulse response is then determined by any 
adapted method. A first step of the method then consists of searching, over said portion 
of the module estimation of the pulse response of the channel, the pulse having a 

25 maximum amplitude. Such an pulse is associated with a reference time tR and is 
considered as corresponding to the main path taken by the signal received by the 
demodulator. In the case where the entire module estimation of the pulse response is 
used, the main path is determined by any adapted method. 

Filtering function FE is then displaced with respect to time tR from an initial time 

30 tl, smaller than tR, to a final time t2, greater than tl and smaller than tR, so that the step 
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of index mi is located at initial time tl before main path 60, and at final time t2 after 
main path 60. 

For a determined position of filtering function FE, the energy of the estimation of 
the pulse response with respect to filtering function FE is obtained by multiplying the 
5 amplitude of each echo by the amplitude of the step containing said echo, by multiplying 
the amplitude of the main path by the amplitude of the step containing the main path and 
by adding the amplitudes thus weighted. In practice, the module estimation of the pulse 
response is stored in the form of a table containing indexes representative of successive 
times, a value of the module estimation of the pulse response being associated with each 

10 index. An analog table representative of filtering function FE is then determined at the 
initial position at time tl. In the case where a single portion of the pulse response is used 
to determine the energy, a portion only of the table representative of the module 
estimation of the pulse response is kept. The energy of the module estimation of the 
pulse response with respect to filtering function FE at initial position tl is obtained by 

15 multiplying, for each index, the values stored in the tables representative of the filtering 
function and of the module estimation of the pulse response, and by adding the obtained 
products. Shifting the filtering function with respect to the module estimation of the 
pulse response then amounts to shifting the indexes of the table representative of the 
filtering function. 

20 The applicant has shown that the use of a filtering function FE such as defined by 

the previously-mentioned characteristics advantageously enables taking into account the 
echoes outside of the guard interval for the positioning of the FFT analysis window while 
avoiding an instability of the method. 

In Figs. 10A, 10B to 13 A, 13B, which illustrate the advantages of the present 

25 invention, filtering function FE is shown, for simplicity, with a smaller number of steps 
than the filtering function of Fig. 7. Each Fig. 10A, 10B to 13 A, 13B shows a portion of 
the module estimation of the pulse response of the channel, and more specifically main 
path 60 and an additional echo (Figs. 10A, 10B to 12 A, 12B) or two additional echoes 
(Figs. 13A, 13B). For each pair of Figs. 10A, 10B, 12A, 12B, and 13A, 13B, the relative 

30 positions between the main path and the additional echo(es) are identical. Figs. 10A, 
10B (respectively 12 A, 12B, and 13 A, 13B) show two positions of filtering function FE 
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with respect to the module estimation of the pulse response and illustrate which position 
of filtering function FE corresponds to a maximum energy. 

Figs. 10A, 10B illustrate the case in which an echo 62 substantially of same 
amplitude as main path 60 is separated therefrom by a duration smaller than the duration 
of the step of index mi, that is, of a duration smaller than the duration of guard interval 
Tg. In this case, whatever the position of filtering function FE between the position 
shown in Fig. 10A and the position shown in Fig. 10B, the energy of the module 
estimation of the pulse response is identical. A conventional result obtained with a 
window having a width equal to guard interval Tg is then found. 

Figs. 11 A, 11B illustrate the case in which an echo 64 is substantially of same 
amplitude as main path 60 is present outside of guard interval Tg, that is, the duration 
between main path 60 and echo 64 is greater than the duration of guard interval Tg. In 
this case, the obtained energy is greater when filtering function FE occupies the position 
shown in Fig. 1 1 A with respect to he position shown in Fig. 1 IB. The positioning of the 
FFT analysis window is thus determined from the position of the filtering function 
corresponding to Fig. 11 A. The present invention thus enables taking into account 
echoes of strong amplitude outside of the guard interval. 

Figs. 12 A, 12B illustrate the case in which an echo 66 having an amplitude 
smaller than the main path is present outside of the guard interval. In this case, the 
obtained energy is greater when filtering function FE occupies the position shown in Fig. 
12A with respect to the position shown in Fig. 12B. This means that the participation of 
an echo 66 of small amplitude outside of guard interval Tg is smaller than the 
participation of main path 60 for the positioning of the FFT analysis window. 

Figs. 13 A, 13B illustrate the case in which two echoes 68, 70 of small identical 
amplitudes are arranged symmetrically on either side of main path 60, the duration 
separating each echo 68, 70 and main path 60 being greater than half of the duration of 
guard interval Tg. In this case, the obtained energy is maximum when filtering function 
FE occupies the position shown in Fig. 13A with respect to the position shown in Fig. 
13B. This means that echoes of same amplitude symmetrical with respect to the main 
path outside of the guard interval have an identical participation for the positioning of the 
FFT analysis window. 
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The present invention thus enables taking into account echoes outside of the 
guard interval for the positioning of the FFT analysis window, while weighting the 
participation assigned by such echoes according to their amplitude and to the deviation 
with respect to the main path, enabling avoiding an instability of the method for 
positioning the FFT analysis window. 

The position of maximum energy obtained from filtering function FE can be used 
to refine the positioning of the interpolation filter implemented to determine the 
estimation of the frequency response of the channel for carriers other than the continuous 
and scattered pilots. 

Of course, the present invention is likely to have various, alterations, 
improvements, and modifications which will readily occur to those skilled in the art. In 
particular, in the example of demodulator of Fig. 6, all units may be modified or replaced 
with appropriate elements. For example, input E of the circuit can directly receive a 
signal centered on approximately 4.5 MHz. The analog-to-digital converter may be 
external to the demodulator. 

The present invention has mainly been described in the context of the digital 
television radio transmission, defined by standard ETSI EN 300 744, VI. 4.1. However, 
the present invention is neither limited to this standard, nor to this field, and may be 
applied in and to any device comprising a COFDM demodulator, be it a television 
receiver or not. For example, the demodulator according to the present invention may be 
used in a portable phone. 

Further, filtering function FE has been described in the context of a digital signal 
processing method. Clearly, the present invention may be implemented in the context of 
an analog signal processing method, the filtering function used being a continuous 
function then corresponding to the previously-described stepped function. 

What is claimed is: 



